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Abstract 

We present the simulation of the time evolution of the distance matrix. 
The result is the node-node distance distribution for various kinds of net- 
works. For the exponential trees, analytical formulas are derived for the 
moments of the distance distribution. 



1 Introduction 

A graph is defined as a set of nodes (vertices) and a set of links between nodes 
(edges) P |2 El EH • By graph evolution or growth we mean subsequent attach- 
ing of new nodes with m edges to previously existing nodes 0. Such growing 
graphs may reflect some features of real evolving networks, e.g. a network of 
collaborators, a network of citations of scientific papers, some biological net- 
works (food chains or sexual relations) or Internet and world-wide-web pages 
with links between them 001311111]. 

The distance between nodes is the shortest number of edges which leads from 
one node to the other. The node-node distance (NND) distribution depends on 
how subsequent nodes are attached. If each node is connected with only one 
of preexisting nodes (m = 1) a tree appears. When a new node is attached to 
several different nodes with m > 1 edges, the growing structure is called a simple 
graph. We may choose nodes to which new nodes are attached in preferential 
way or randomly. In the latter case we deal with exponential networks. If the 
probability of choosing a node is proportional to its degree (e.g. to the number 
of its nearest neighbors) the growing structure is called scale-free or Albert- 
Barabasi networks |10j . 

In this paper, the numerical algorithm for the network growth — basing on 
distance matrix evolution — is presented both, for exponential and scale-free 
networks (m — 1,2) [111 112) . The NND distribution and its characteristics are 
calculated. For the exponential trees the iterative formulas for n-th ordinary 
moments arc derived. 



2 Computer simulations 

A graph with edges of unit length may be fully characterized by its distance 
matrix S, an element Sij of which is equal to the shortest path between nodes 
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i and j. This matrix representation is also particularly useful when computer 
simulations for graph evolution are applied. 

Attaching subsequent node with one edge (to = 1) to previously existing 
network of N nodes corresponds to adding a new (N + l)-th row and a new 
column to N x N large distance matrix S. The distance from newly added 
(N + l)-th node to all others via selected — labeled as p — node is larger 
by one than the distance between p-th and all others. Thus, new (N + l)-th 
row/column is a simple copy of p-th row/column but with all of its elements 
incremented 

V 1 < i < N : s N+1 (N + l,i) = s N+1 (i, N + l) = s N {p, i) + 1. (1) 

Similarly, when new node is attached to the network with two edges (to = 2) 
to two different nodes — which are labeled as p and q — but the distance from 
all other nodes i to the newly added (N + l)-th is one plus the smaller distance 
between p-i or q-i nodes pairs ^2] : 

VI < i < N : s N+1 (N+l,i) = s N+1 (i,N+l) = mm(s N (p,i),s N (q,i)) + l. (2) 

In the case mentioned above of the growth of the simple graphs also the 
reevaluation of distances between nodes i and j must be done to check if adding 
a new node provides the shortcut |12j : 

Vl<i,j <N:s N+1 (i,j) = 
min (sN(i,j),SN(i,p) + 2 + s N (q,j), s N (i,q) + 2 + s N (p,j)). 

In both cases diagonal elements of new row/column are zero |11M12) : 

s N +i(N + 1,JV + 1)=0. (4) 

Selecting rows/columns (nodes to which we attempt to add a new node) may 
be random or preferential. In the latter case an additional evolving vector is in- 
troduced, which contains the node labels. These labels occur as vector elements 
with a probability proportional to the degree of the node. Random selection of 
elements of such a vector correspond to Albert-Barabasi construction rule. The 
procedure is known as the Kertesz algorithm |13| . 



3 Analytical calculations 

Let us define n-th moments of the NND distribution for all distances 

N N 

r N = [{s n (i,j)}] ^^^[^(m)], (5) 

i=i j=i 

and only for non-zero distances 

JV N 

d n N = [(s"(i,j))] = (6) 
\ I i=l j=\ 

where {••■}, (•••) and [• • •] denote an average over N 2 matrix elements, an 
average over N(N — 1) non-diagonal matrix elements, and an average over N lun 
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Table 1: The mean distance d(N) — a In N + b for different evolving networks. 





exponential 


exponential 


scale- free 


scale-free 


m 


1 


2 


f 


2 


a 


2.00 


0.672 


1.00 


0.462 


b 


-2.84 


0.296 


-0.08 


0.889 



Table 2: The dispersion a 2 (N) = c In N+d for the exponential and the scale- free 
trees (m = 1). 



exponential 


scale-free 


c 2.00 


1.00 


d -1.44 


-1.64 



independent realizations of the evolution process (matrices), respectively. Mo- 
ments iJUJ and Jfjjl for n = 1 are sometimes called the network diameter. Both 
double sums in r.h.s. of Eqs. © and © are equal, due to obvious fact, that 
s(i,i) = 0. That allows to derive simple dependence between averages {• ■ • } 
and (■■•): 

Nil = (N- l)d n N . (7) 

For the exponential trees — basing on s(i, i) = and distance matrix symmetry 
s ihj) = — we are a ble to construct iterative equations for f^.j as 

dependent on £% (k = 1, • • • ,n) 

N+1N+1 N 

(N + lft n N+1 = £ E t 5 "^)] - n2£ n + 2 E i 1 + M (8) 

i— 1 J — 1 z— 1 

where q is the number of the randomly selected row/column of the distance 
matrix S. Combining Eq. J3J) with Eq. J7J) gives the desired iterative formula: 

_ {N + 2){N-l) 2 2(JV-1) y^M fe 

dN+1 ~ (N + l)N dN + WTT + (N + 1)N^ [kj 



4 Results and conclusions 

For the trees the mean of the NND djy and its dispersion a 2 = d 2 N — {d]^) 2 
grow logarithmically with N (see Tabs. ^ 13 HQ- F° r the graphs only the 
first cumulant (the average of the NND d} N ) grows logarithmically (see Tab. |2J| 
|12j . Such a slow increase of d]^ with number of network nodes is known as the 
small- world effect |T3] . 



The histogram of NND is presented in Fig. ^ As we expected the NND 
for the graphs are more condensed than the NND for the trees, as well as, the 
scale- free graphs (trees) are more condensed than the exponential graphs (trees). 

Knowing the moments d^ — the averages of the n-th powers of the non- 
diagonal distance matrix elements © — allows to build all statistical parame- 
ters which characterize the NND distribution, e.g. the average distance d, the 
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N=1 000, N run =1 3 (m=1 ), N mn =1 <f (m=2) 
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Figure 1: The NND distribution for different types trees and graphs. 



exponential trees, N run =10 
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Figure 2: The main moments d 1 ^ (k = 1, ••• ,4) for the exponential trees given 
by Eq. (lines) and from the direct simulations (symbols). The latter are 
averaged over N inn = 10 4 independent evolution process realizations. 
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node-node distance distribution characteristics for exponential tress 
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Figure 3: The NND distribution characteristics for the exponential trees as 
derived from iterative Eq. (jty. 



distance dispersion a 2 , its skewness 

_ i N — u.i x a N -r -\'< \ 

and kurtosis 



d% - M^dj, + 2(d 1 N ) 3 



_ d%- 4d|4 + M 2 ^) 2 - 3(d] v ) 4 
K4 = -3 . 

The values of such characteristics of the NND for the exponential trees obtained 
via Eq. © are presented in Figs. [2land[3] For trees, the distributions are similar 
to the Poisson distribution (see Fig. QJ. However, even for large N the skewness 
and kurtosis do not vanish as one may expect for the normal distribution |15|. 
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